Nanofibers with Bragg gratings from equidistant holes 
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We study nanofibers with Bragg gratings from equidistant holes. We calculate 
analytically and numerically the reflection and transmission coefficients for a 
single grating and also for a cavity formed by two gratings. We show that the 
reflection and transmission coefficients of the gratings substantially depend on 
the number of holes, the hole length, the hole depth, the grating period, and 
the light wavelength. We find that the reflection and transmission coefficients 
of the gratings depend on the orientation of the polarization vector of light 
with respect to the holes. Such a dependence is a result of the fact that the 
cross section of the gratings is not cylindrically symmetric. 
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1. Introduction 

A fiber Bragg grating (FBG) is a structure with a peri- 
odic variation in the refractive index in a small section 
of an optical fiber. A FBG can be formed from multi- 
ple layers of alternating materials with varying refractive 
index, or by a periodic variation of some characteristic 
of the fiber core. The change in the refractive index in 
the axial direction of the fiber causes a partial reflec- 
tion of an optical wave. For waves whose wavelength is 
close to twice the grating period, the many reflections 
combine with constructive interference, and the struc- 
ture acts as a high-quality reflector. The first FBG was 
demonstrated by Hill et al. using a visible laser propa- 
gating along the fiber core [1] . Meltz et al. demonstrated 
the much more flexible transverse holographic technique 
where the laser illumination came from the side of the 
fiber [2]. This technique uses the interference pattern of 
ultraviolet laser light to create the periodic structure of 
the Bragg grating. FBGs have been widely used in many 
applications [3-7]. The primary application of FBGs is 
in optical communications systems. They are specifically 
used as notch filters. They are also used in optical mul- 
tiplexers, demultiplexers, and optical fiber sensors. 

Due to recent developments in the taper fiber technol- 
ogy, thin fibers can be produced with diameters smaller 
than the wavelength of light [8, 9]. The essence of the 
technology is to heat and pull a single-mode optical fiber 
to a very small thickness, maintaining the taper condi- 
tion to keep adiabatically the single-mode condition [10]. 
Due to tapering, the original core is almost vanishing. 
Therefore, the refractive indices that determine the guid- 
ing properties of the tapered fiber are the refractive in- 
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dex of the original silica clad and the refractive index of 
the surrounding vacuum. Such subwavelength-diameter 
vacuum-clad silica-core fibers arc called nanofibers. Sev- 
eral methods for trapping and guiding neutral atoms out- 
side a nanofiber have been proposed and studied [11- 
15]. A tapered fiber with an intense evanescent field 
can be used as an atomic mirror [16]. The genera- 
tion of light with a supercontinuum spectrum in thin 
tapered fibers has been demonstrated [9]. The evanes- 
cent waves from zero-mode metal-clad subwavelength- 
diameter waveguides have been used for optical obser- 
vations of single- molecule dynamics [17]. Efficient chan- 
neling of emission from a few atoms into guided modes 
has been realized [18]. The correlations between pho- 
tons emitted from a few atoms have been measured [19]. 
Thin fiber structures can be used as building blocks in 
future atom and photonic micro- and nanodevices. Thin 
fibers can also be used as efficient nanoprobes for atoms, 
molecules, and quantum dots. 

Recently, it has been proposed to combine the FBG 
technique with the nanofiber technique to obtain a hy- 
brid system [20]. In such a system, two FBG mirrors 
form a cavity transmitting and reflecting the guided field 
of the nanofiber. Therefore, the interaction between the 
field and the atoms in the vicinity of the fiber is enhanced 
not only by the transverse confinement of the field in the 
fiber cross-section plane but also by the longitudinal con- 
finement of the field in the space between the mirrors. An 
advantage of the nanofiber cavity is that the field in the 
guided modes can be confined to a small cross-section 
area whose size is comparable to the light wavelength. 
Another advantage of the nanofiber cavity is that the cav- 
ity guided field can be transmitted over long distances for 
communication purposes. The nanofiber cavity can find 
applications in the merged areas of fiber optics, cavity 
quantum electrodynamics, ultracold neutral atoms, and 
electromagnetically induced transparency [21-26]. It has 
been shown that a nanofiber cavity with a large length 
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(on the order of 10 cm) and a moderate finesse (about 
30) can significantly enhance the group delay of a guided 
probe field [20] and substantially enhance the channeling 
of emission from an atom into the nanostructure [27]. 
The effect of an atom on a quantum guided field in a 
weakly driven FBG cavity has been studied [28]. 

In order to produce a fiber cavity, one may splice com- 
mercial FBGs to a fiber. However, this technique is not 
suitable for making a nanofiber cavity because it creates 
substantial losses in the tapering region. Meanwhile, the 
UV irradiation technique for making FBGs uses the pho- 
tosensitivity of the Ge/Ge0 2 -doped core of optical fibers. 
However, this technique does not work for nanofibers be- 
cause the original core is almost vanishing while the origi- 
nal clad acts as the core of the nanofiber. Recently, a new 
technique for making nanofiber cavities has been demon- 
strated [29]. This technique uses focused ion beams [30] 
to drill periodic holes on a nanofiber. The typical fi- 
nesse of a nanofiber cavity produced by the focused-ion- 
beam milling technique is F ~ 20-120 and the typical 
on-resonance transmission is ~ 30-80%. Such nanofiber 
cavities can find applications in the merged areas of fiber 
optics and cavity quantum electrodynamics [20-26]. 

In this article, we study nanofibers with Bragg gratings 
from equidistant holes. We calculate analytically and 
numerically the reflection and transmission coefficients 
for a single grating and also for a cavity formed by two 
gratings. We find that the reflection and transmission 
coefficients of the gratings depend on the orientation of 
the polarization vector of light with respect to the holes. 

The article is organized as follows. In Sec. 2 we de- 
scribe the model and present the basic equations of the 
coupled-mode theory for the system. In Sec. 3 we calcu- 
late the reflection and transmission coefficients of a single 
grating. In Sec. 4 we calculate the reflection and trans- 
mission coefficients of a cavity formed by two gratings. 
Our conclusions are given in Sec. 5. 

2. Model and basic equations 

We consider a nanofiber with a grating of equidistant 
pairs of lateral identical holes (see Fig. 1). The nanofiber 
has a cylindrical silica core of radius a and of refractive 
index n\ — 1.45 and an infinite vacuum clad of refractive 
index n-i — 1. In view of the very low losses of silica 
in the wavelength range of interest, we neglect material 
absorption. We use the Cartesian coordinates {x, y, z) 
and the cylindrical coordinates (r, tp, z), with z being the 
axis of the fiber. 

The holes are carved on the opposite lateral sides of 
the fiber [see Fig. 1(a)]. They are paired, one left and 
one right. Such a grating can be produced [29] by the 
focused ion beam milling method [30]. The cross section 
of the nanofiber in the absence of holes is a disk [see 
Fig. 1(b)]. The medium in this area is silica, with the 
refractive index n\ — 1.45. The cross section of a hole 
has the form of a circular segment [see Fig. 1(c)]. The 
width d of the circular segment is the depth of the hole. 



Each hole is symmetric about the axis x. The holes of 
every pair are symmetric to each other with respect to 
the axis y. The symmetry axes x and y and are called 
the principal axes. The length of a hole along the fiber 
axial direction z is denoted by h. The spatial period of 
the grating of holes is denoted by Ag ■ Since the holes are 
simply vacuum voids, the refractive index of the medium 
in the holes is ni = 1. 
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nanofiber cross section of lateral holes 



Fig. 1. (a) Nanofiber with a grating of equidistant pairs of 
lateral identical holes, (b) Cross section of the nanofiber 
in the transverse plane, (c) Cross section of a pair of 
holes carved into the fiber on the opposite lateral sides. 

To describe the field, we use the basis consisting of 
the modes of the nanofiber in the absence of the holes. 
In an unperturbed (ideal) fiber, the electromagnetic field 
can be decomposed into guided and radiation modes [31]. 
The guided modes have the evanescent behavior on the 
outside of the core. They can travel in the waveguide 
without loss of power, provided that losses in the di- 
electric material are ignored. Meanwhile, the radiation 
modes are oscillatory at large distances from the fiber 
and do not have the evanescent behavior. They cannot 
be normalized to a finite amount of power. 

Regarding the guided modes, we assume that the 
single-mode condition [31] is satisfied for a finite band- 
width around a central frequency of interest. We label 
the fundamental guided modes HEn in this bandwidth 
by the index fj, = (ui,f,l), where uj is the frequency of 
the mode, / = +,— denotes the forward or backward 
propagation direction and I — +, — stands for the coun- 
terclockwise or clockwise rotation of polarization. The 
propagation constant j3 — /3(u>) for a guided mode is 
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determined by the fiber eigenvalue equation [31]. The 
profile functions £^ = e^e* and = h^e ilv of 
the electric and magnetic components, respectively, of 
the field in the guided mode [i = (u>, f, I) are given in 
Appendix A. 

Regarding the radiation modes, the longitudinal prop- 
agation constant f3 for a given frequency ui can vary con- 
tinuously, from — k to k, with k — ujjc being the wave 
number. We label the radiation modes by the index 
v = (uj,(3,m,l), where m = 0, ±1,±2, ... is the mode 
order and I = +, — is the mode polarization. The profile 
functions £ v = e^e im f and U v = h^e im ^ of the elec- 
tric and magnetic components, respectively, of the field in 
the radiation mode v — (oj, (3, m, I) are given in [27, 31]. 

We assume that the field has a single frequency lj. 
The orthogonality of the modes of the unperturbed fiber 
is expressed by the formula 

X - [Re[u z - (Si xn*,)]dr = P i 6 ii >. (1) 

Here u z is the unit vector for the direction z, and 5u> is 
the Kronecker delta symbol for discrete indices i and i' 
of guided modes, the Dirac delta function for continuous 
indices i and i' of radiation modes, and zero if one of the 
indices i and i' is discrete and the other one is continuous. 
We have used the notation J dr = J Q x dp J °° r dr for the 
integral over the fiber transverse plane. The notation Pi 
stands for the power of the mode i. For convenience, we 
normalize the mode functions to the same power P, that 
is, we choose Pi = P. 

We now consider the field in the presence of the holes. 
It is our purpose to expand the field in terms of the modes 
of the ideal fiber. We separate the field vectors into the 
transverse and longitudinal parts designated by the sub- 
scripts t and z as 

E = E t + E z , 

H = H i + H z . (2) 

Due to the completeness of the set of transverse mode 
functions, we can expand the transverse components E t 
and H t of the electric and magnetic parts of the field in 
terms of the ideal modes as 

a 

= £>+ - a-)H at . (3) 

a 

Here a+ and a~ are the amplitudes of the field in the for- 
ward and backward modes, respectively, £ a and H a are 
the profile functions for the electric and magnetic parts of 
the forward modes, and a = fi, v is the common label for 
the guided and radiation modes. The longitudinal com- 
ponents E z and H z of the electric and magnetic parts of 



the field can be written as 

a 

H, = 5>^ (4) 

a 

where the functions n = n(x, y, z) = n(r, ip, z) and no — 
n (x,y) = no(r,p) are the refractive-index distributions 
of the perturbed and unperturbed fibers, respectively. 

According to the coupled- mode theory [31], the ampli- 
tudes a+ and a~ are governed by the equations 

^ = ifi a a+ +i^2(U aa ,a+, + V aa >a~,), 

a' 

= -i& a a~ - i^(V aa *a+, + U aa ,a~,), (5) 

a' 

where 

U aa > = ^ J(n 2 - nl)[(£ at ■ S* a , t ) 
+ (nl/n 2 )(£ az -£* a , z )}dr, 

V aa , = ^ J(ni-nl)l(£ at -£* a , t ) 

-(nl/n*)(£ az -£* a , z )]dv (6) 

are the mode-coupling coefficients. We note that the 
coupled-mode equations (5) involve the complete set of 
guided and radiation modes and are exact. 

We now assume that the coupling between the guided 
modes and the radiation modes is negligible, that is, the 
grating of the holes along the fiber is lossless. When we 
omit the radiation modes from Eqs. (5), we obtain the 
following propagation equations for the coupled guided 
modes: 

da + 

= i(3 ia + + iJ2{Uwa+ + 
^ = -ifrai-i^2(V w a+ + U w a^). (7) 

Here I is the guided-mode index and takes the value 1 
or —1, corresponding to the left or right circular polar- 
ization, respectively. We emphasize that the omission 
of the radiation modes from the coupled-mode equations 
is the only approximation used in this paper. This ap- 
proximation is valid only when the losses are not serious, 
that is, when the characteristic size of the holes is small 
as compared to the fiber radius. It is interesting to note 
that, according to [29], the losses are not too serious (20- 
70%) even when the length and depth of the holes are 
substantial fractions of the fiber radius. For such losses, 
our analytical and numerical results shown below are not 
quantitatively valid. However, they can still be used to 
qualitatively understand the underlying physics and to 
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explain many features of the nanofiber cavities produced 
by the focused-ion-beam milling method [29]. 

Note that j3\ — 0-i = ft. Due to the symmetry of 
the mode functions and the hole shape, we have U± 1 = 

17-1,-1 = U, V 1}1 = V-!,-! = V, [/!,_! = U-1,1 = U C , 

and Vi,-i = V-i.i = V c . In addition, we can show that 
U = U*,V =V*,U C = U*, and V c = V*. Hence, Eqs. 
(7) become 



da* 
dz 

da* 1 
dz 
da i 
dz 

daZi 
dz 



i(0 + U)a + + iVaZ + iU c a\ + iV c aZ 1 , 
i{0 + U)aZ 1 + iVaZi + iU c a\ + iV c aZ , 
— i(0 + U)aZ — iVa + — %U c aZ\ — iV c aZ\, 
-i{/3 + U)aZ 1 - iVati - iU c a{ - iV c af.(8) 



We note that the coefficients U c and V c characterize the 
coupling between the two different circular polarizations 
I = ±1. Such cross coupling arises because the cylindrical 
symmetry of the fiber is broken by the presence of the 
holes on the lateral sides. 

We introduce the notations a x = (a\ + a_i)/\/2 and 
a y = ( a i — fl-i)/iV2| which stand for the amplitudes of 
the fields that are linearly polarized along the x and y 
directions. Then, Eqs. (8) split into two different sets of 
coupled equations, namely 



and 



da* 

dz 
da~ 

dz 

da+_ 
dz 

da~ 
dz 



i(p + U x )a+ - 
-i{p + U x )a~ x 



iV x a x , 



iV x a + , 



(9) 



i(/3 + U y )a+ +iV v a~, 
-i(P + U y )a~ -iV y a + . 



(10) 



Here we have introduced the notations 



U x = 

v x = 

Uy = 

v„ = 



u + u c , 
v + v c , 

U-Uc, 

V -v c . 



(11) 



Equations (9) and (10) show that the x and y polariza- 
tions are not coupled to each other. Consequently, the 
x and y polarizations can be considered as the principal 
polarizations of the field in the nanofiber in the presence 
of the holes. Equations (11) show that the principal po- 
larizations x and y have different coupling coefficients. 

3. Single gratings 

In this section, we calculate the reflection and transmis- 
sion coefficients of a single grating. We start our calcula- 
tions by considering a single step of a grating, that is, a 



single pair of holes. We assume that the two holes stretch 
from the point z = to the point z = h along the fiber 
axis. In the interval (0, h), the coupling coefficients U a 
and V a , where a = x,y, are constant. For the boundary 
condition a~(h) = 0, the solutions to Eqs. (9) and (10) 
are found to be 



at{z) 
a~(z) 



A a {K a cos[K a {z-h)} 
+ i(/3 + U a )sm[KAz-h)]}, 
-iA a V a sm[K a {z - h)], 



where 



K„ 



y/W + u.y-v*. 



(12) 



(13) 



The coefficient A„ is determined by the input field am- 
plitude a+(0). 

From Eqs. (12) and (13), the reflection coefficient 
R a = a~(0)/a + (0) and the transmission coefficient T a = 
a + (h)/a + (0) of the single pair of holes for the field with 
the principal polarization a — x,y are found to be 



iV a sin{K a K) 



K a cos(K a h) -i{P + U„) sm(K a h) ' 

T = — (Id) 

K a cos(K a h)-i(l3 + U a )sm(K a h)' v ; 

It is clear that the reflection and transmission coefficients 
of the pair of holes depend on the orientation of the po- 
larization vector of the field. The reason is that, when the 
total cross section of the pair of holes is not cylindrically 
symmetric, the cross-coupling coefficients U c and V c be- 
come nonzero. Therefore, we have U x ^ U y and V x ^ V y 
[see Eqs. (11)], and hence K x ^ K y [sec Eq. (13)]. This 
leads to R x ^ R y and T x =fi T y [see Eqs. (14)]. 

We now consider a grating of TV equidistant identical 
pairs of holes stretching from the position z — to the 
position z = Lq = Nh + (N — l)s. Here s — Xq — h is 
the length of the region that separates two adjacent pairs 
of holes. We consider a principal polarization a = x, y. 
The solution of Eqs. (9) or (10) for the first step of the 
grating (the first pair of holes) can be written as 



where 



a+_{h) 
a a {h) 



M = 



M 



a+(0) 
M0) 



Mn M12 
M 2 i M 22 



(15) 



(16) 



is the transfer matrix for a single step, with the matrix 
elements 



M 22 
Mi 2 
M 2i 



cos(Kcrh) + i — ■ jpr— sin(K a h), 
cos(if CT /i) — i — — s\n(K a h), 



i—sin(K a h), 
-ijjr- sm(K a h). 



(17) 
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Note that M 12 = -M 21 and M X1 M 22 - M 12 M 21 = 1. 

In the region from z — h to z = h + s = Xq, which sep- 
arates the first and second steps, the fiber is unperturbed 
and, hence, the field propagates as an ideal guided-mode 
field with the propagation constant /3. Therefore, we have 



a+_(h 
a a (h 



a+Jh) 
a a (h) 



where 







-i/3s 



(18) 



(19) 



is the propagator for the ideal guided modes of the fiber. 

When we apply the above procedures to the N steps 
of the grating and combine the results, we find the input- 
output relation 



a+(L)\ _ wW /o+(0) 
a~(L) J- W lo-(0) 



where 



W (JV) = M(FM) 



N-l 



(20) 



(21) 



is the total transfer matrix for the grating. We note 
that detWW 



W- 



(N) 



= 1 and 



In addition, we have W^- 1 ) = M. 



The explicit expressions for the elements of the total 
transfer matrix are found to be 



sinh(A^) 
M\i — . , „ e 



II 



(N) 



22 



II 



(AT) 



12 



Mi 



Mai 



sinh 9 a 
smh(N6 a ) 
1 smh9„ 
sinh(N6 a ) 

sinh 9 a 
sinh(iV0 CT ) 

sinhfl CT : 



_ iga sinh[(JV-l)fl g ] 

sinh 9 a 
ifjs smh[(N - l)9 a ] 
sinh No- 



where 



= ln( J D CT + /Df^T), 



(22) 



(23) 



with D a = |(M n e^ s + Af 22 e~ l/3s ), that is, 



D a = cos(K a h) cos(/3s) — sm.{K a h) sin(/3s). 

Note that 9 a is, in general, a complex number. 

The reflection and transmission coefficients of the grat- 
ing are given by R a N) = W$P /W$P = /W^ 

and T„ — 1/W 22 \ respectively. The explicit expres- 
sions for these coefficients are found to be 



T (N) 



R a sinh(iV6' (T ) 



smh(N9 a ) - T a e lfis sinh[(7V - l)9 a ] 

T a sinh 9 a 
sinh(7V6> CT ) - T a e^ s sinh[(7V - l)9 a ] ' 



(25) 



Here R a 



R 



(N=l) 



T, 



(N=l) 



M 12 /M 2 



-M 21 /M 22 and 



1/M 22 are the reflection and transmis- 
sion coefficients for a single step and arc given by Eqs. 
(14). Equations (25) show that, like the coefficients R a 

and T a , the coefficients Ra and T a N ^ depend on the 
principal polarization a — x, y of the field. In addition, 

since 9 a is, in general, a complex number, Ra ^ and T a N ^ 
may oscillate when a parameter, such as the number of 
holes N, the grating period Ag, the hole depth d, the 
hole length h, and the wavelength of light A, varies. 

We note that, when we use Eq. (21), we can derive the 
recurrence formulae 

l( - N}2 R e 2%l3s 



R (N+1) = R (N) + 



1 - E^R^e^P' 

i N )rp AB S 

7 ± 



(26) 



We can also obtain these formulae in the framework of ray 
optics. The explicit expressions (25) for Ra and T a N ^ 
are in agreement with the recurrence formulae (26). 

We calculate the reflection and transmission coeffi- 
cients Ra and T a N ^ from expressions (25). We plot 

in Figs. 2-6 the reflectivity \Ra | 2 as a function of the 
number of steps N , the grating period Ag, the hole depth 
d, the hole length h, and the wavelength of light A. Com- 
parison between the solid curves (for the x polarization) 
and the dashed curves (for the y polarization) shows that 
the different principal x and y polarizations of light lead 
to different reflectivities of the grating. 



(a) 



0J 
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(c) 



300 



Number of steps N 



Fig. 2. Reflectivity | 2 as a function of the number 

of steps N. The solid and dashed curves are for the prin- 
cipal x and y polarizations, respectively. The wavelength 
of light is A = 852 nm. The grating period is Ag = 363 
nm (a), 364.5 nm (b), and 366 nm (c). The fiber radius is 
a = 290 nm. The hole length and hole depth are h = 150 
nm and d = 100 nm, respectively. 

We observe from Fig. 2 that the dependence of \Ra N ^\ 2 
on N can be either monotonic or oscillatory depending 



on the parameters of the fiber, the grating, and the light 
field. The reason for the oscillatory behavior is that the 
parameter 9 a is, in general, a complex number, and there- 
fore the hyperbolic functions sinh(A r 6' cr ) and cosh(A r 6' cr ) 
in expressions (25) may become oscillatory functions. In 
the case of Fig. 2(a), Ifi^'l 2 (dashed line) dominates 
l-R^I 2 (solid line). In the case of Fig. 2(c), |iZjfe*°| a 
(solid line) dominates l-Ry^l 2 (dashed line) when N is 
large enough. In the case of Fig. 2(b), both |i?£ JV ' ) | 2 (solid 

line) and l-R^I 2 (dashed line) achieve similar significant 
values when N is around 200. 
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Fig. 3. Reflectivity \R^\ 2 as a function of the grating 
period Aq. The solid and dashed curves are for the prin- 
cipal x and y polarizations, respectively. The number of 
steps is N = 60 (a), 120 (b), and 240 (c). Other param- 
eters are as in Fig. 2. 
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Fig. 4. Reflectivity \R^ | 2 as a function of the hole depth 
d. The solid and dashed curves are for the principal x 
and y polarizations, respectively. The number of steps is 
N = 60 (a), 120 (b), and 240 (c). The grating period is 
Xg =364.5 nm. Other parameters are as in Fig. 2. 
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Fig. 5. Reflectivity lii^! 2 as a function of the hole 
length h. The solid and dashed curves are for the prin- 
cipal x and y polarizations, respectively. The number of 
steps is N = 60 (a), 120 (b), and 240 (c). The grat- 
ing period is Xq =364.5 nm. Other parameters are as in 
Fig. 2. 
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Fig. 6. Reflectivity l-R^I 2 as a function of the wave- 
length of light A. The solid and dashed curves are for the 
principal x and y polarizations, respectively. The num- 
ber of steps is TV = 60 (a), 120 (b), and 240 (c). The 
grating period is Xg =364.5 nm. Other parameters are 
as in Fig. 2. 



Figures 3-6 show that the reflectivity | 2 of the 

grating oscillates with increasing Xq, d, h, or A. We 
observe that l-R^I 2 can achieve significant values when 
the parameters Xg, d, h, or A are appropriate. Figures 3— 
6 show that the region where \R X N | 2 achieves significant 
values is different from the region where l-R^I 2 does so. 
The magnitude and width of the maximum of \Ry | 2 are, 
in general, larger than those of the maximum of |i?x^| 2 - 
Figures 3(c), 4(c), 5(c), and 6(c) show that, for N = 240, 



both \R. 



W|2 



and lit' 



W|2 



can achieve similar significant 
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values when Xg = 364.5 nm, d = 100 nm, h — 150 nm, 
and A = 852 nm. 

We recognize that the length h and the depth d of 
the holes used in the above numerical results are sub- 
stantial fractions of the fiber radius. For such a grating, 
the losses to radiation modes are not small. Therefore, 
our analytical and numerical results are not quantita- 
tively valid. However, they can still be used to qualita- 
tively understand the underlying physics and to explain 
many features of the nanofiber cavities produced by the 
focused- ion-beam milling method [29]. 

4. Nanofiber cavity 

We consider a cavity formed by a nanofiber with two 
Bragg gratings (see Fig. 7). We label the gratings by the 
index j, where j = 1,2. Each of the gratings can be ei- 
ther a single grating or a set of multiple gratings. In the 
case of a multiple grating, the reflection coefficients for 
the left and right sides of the grating may differ from each 
other. Therefore, to be general, we denote the reflection 
coefficient for a light field with a principal polarization 
a = x,y, incident onto the left or right side of grating 
j = 1,2 by the notation R a j or R' a j, respectively. The 
transmission coefficient for grating j is T a j and does not 
depend on the incidence side of the grating. The distance 
between the two gratings of the cavity is denoted by L. 
When we use ray optics and combine multiply reflected 
beams, we can easily calculate the reflection and trans- 
mission coefficients of the set of the two gratings, that is, 
for the whole nanofiber cavity. 



grating 1 



grating 2 





Fig. 7. Cavity formed by a nanofiber with two Bragg 
gratings. 

The reflection coefficients R a and R'^ for the left and 
right sides, respectively, of the set of the two gratings are 
found to be 



RL = R'„ 



cr2 



T 2 al R a2 e^ L 
l-R'^R^e 2 ^' 

T 2 D/ 2 ij3L 
1 cr 2 n al e 

l-R' al R a2 e^- 



(27) 



Meanwhile, the transmission coefficient T„ for the set of 
the two gratings is found to be 



T„ = 



T al T a2 e^ L 
l-R' al R a2 e 2 ^ L ' 



(28) 



Hence, the transmission of the cavity formed by the two 
gratings is 



\Ta- 



?\T« 2 \ 



(1 - Ki||i? ff2 |) 2 + 4| J R al ||.R (72 |sur 2 6 CT ' 



(29) 



where 6 CT = (<j>' al + </> CT 2)/2 + (3L, with <p' al and <i> a2 being 
the phases of the complex reflection coefficients R' al and 
R a2 , respectively. 

We consider the case where the separation between the 
two gratings of the cavity is much larger than the light 
wavelength, that is, L ^ A. We assume that the fre- 
quency <jj of the light field is in a small interval around a 
central frequency u>q. In this case, we can use the expan- 
sion (3(uj) = po + 5/vg, where /3o = /3(wo) is the central 
value of the fiber propagation constant, 5 — u> — luq is the 
detuning of the field, and v g — (dP/du))~ 1 \ LU=UJo is the 
group velocity of the guided light field. Then, we have 



6 CT = (0^ + CT2 )/2 + p L + LS/v g 



(30) 



When the separation L between the two gratings is suffi- 
ciently large as compared to their thicknesses, the mod- 
uli and phases of the complex reflection and transmis- 
sion coefficients of the gratings vary very little with 
the frequency in the range of a few times of nv g /L 
around u>q. Then, we find from Eqs. (29) and (30) 
that the separation between the resonances is Aw rcs = 
TTVg/L and the width of the resonances is 5u> rcs = 
(1 - \Rcri\\Ra 2 \)v g /{^\Rcyi\\R a2 \L). Hence, the finesse 
of the cavity is estimated to be F a = Aw res /&j ros = 
7r-^/ l^o-i 1 1 1/(1 — |-Rcti||-R<t2|)- It is clear that the cav- 
ity transmission \T a \ 2 and the cavity finesse F„ depend 
on the principal polarization a = x, y of the field. 
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Fig. 8. Cavity transmission \T a \ 2 as a function of the de- 
tuning 5. The solid and dashed curves are for the princi- 
pal x and y polarizations, respectively. The cavity length 
is L = 7.6 mm. The two FBG mirrors are identical. The 
number of steps in each mirror is iV = 60 (a), 120 (b), 
and 240 (c). The grating period is Xg =364.5 nm. Other 
parameters are as in Fig. 2. 
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We plot in Fig. 8 the cavity transmission |T CT | 2 as a 
function of the detuning S. Comparison between the solid 
curves (for the x polarization) and the dashed curves (for 
the y polarization) show that the positions and widths of 
the resonances are different for the different principal x 
and y polarizations of the light field. 



0.8-, 




Detuning 5 [GHz] 

Fig. 9. Tuning profile of the cavity transmission \T\ 2 = 
\T X \ 2 cos 2 tp + \T y \ 2 sin 2 tp for tp = n/4. Other parameters 
are as for Fig. 8. 

We emphasize that Eqs. (27)-(29) are valid only when 
the polarization of the light field is linear and is oriented 
along the principal symmetry axis x or y. When the 
polarization of the light field is linear and is oriented 
at an angle tp in the xy plane, the input field can be 
presented in the form Ej n = A(£ x cos tp+£ y simp), where 
£ x = (£\ + £-i)/s/2 and £ y = (£% — £_i)/iv2 are the 
profile functions of the guided modes with the principal 
x and y polarizations, respectively, and A characterizes 
the field strength. Then, the reflected and transmitted 
fields are found to be 

Erf = A(£ X R X cosp + £ y R y sintp), 

Etr = A(£ X T X cos tp + £ y T y sin p) . (31) 

Equations (31) show that, when R x ^ R y and T x ^= T y , 
the directions of the polarization vectors of the input, 
reflected, and transmitted fields are, in general, different 
from each other. Consequently, the polarization vector 
of the field may be rotated after being reflected from or 
transmitted through the cavity. Such a rotation of the 
polarization vector of the field is a result of the breaking 
of the cylindrical symmetry of the fiber. The rotation 
vanishes when ip = or 7r/2, that is, when the polariza- 
tion vector of the light field is oriented along the principal 
axis x or y. 

Due to the mode orthogonality, the power of the field is 
the sum of the powers of the individual x- and y-polarized 
modes. Consequently, the reflectivity and transmission of 
the field in terms of the power are given by 

|i?| 2 = \R x \ 2 cos 2 p+ \R y \ 2 sin 2 tp, 

\T\ 2 = \T X \ 2 cos 2 p + \T y \ 2 sin 2 tp. (32) 



In general, the tuning profile of the cavity transmission 
|T| 2 contains the features of both the profiles of \T X \ 2 and 
\T y \ 2 . We note that the discussions around Eqs. (31) and 
(32) are valid not only for the case of a cavity but also 
for the case of a single grating. 

We plot in Fig. 9 the tuning profile of the cavity 
transmission \T\ 2 — \T X \ 2 cos 2 ip + \T y \ 2 sin 2 ip for the an- 
gle tp = 7r/4 of the polarization orientation of the light 
field. All other parameters are the same as those for Fig. 
8. The figure shows clearly that the profile of |T| 2 has 
double-peak structure. In such structure, one peak cor- 
responds to the resonance in the profile of \T X \ 2 and the 
other peak corresponds to the resonance in the profile of 
\T y \ 2 . Comparison between parts (a), (b), and (c) of the 
figure shows that the tuning profile of the transmission 
of the cavity depends on the number of holes. 
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Fig. 10. Tuning profile of the cavity transmission |T| 2 = 
| T x | 2 cos 2 tp + |T y | 2 sin 2 (p for different values tp = (a), 
7r/6 (b), 7r/4 (c), 7r/3 (d), and it/2 (e) of the orientation 
angle of the polarization vector of the light field. The 
number of steps is N — 210. Other parameters are as for 
Fig. 8. 

We plot in Fig. 10 the tuning profile of the cavity 
transmission |T| 2 = \T X \ 2 cos 2 p + \T y \ 2 sin 2 tp for differ- 
ent values of the orientation angle tp of the polarization 
vector of the light field in the xy plane. It is clear from 
the figure that the tuning profile of \T\ 2 substantially de- 
pends on the angle p. The profile has a single peak when 
tp = (a) or 7r/2 (e), that is, when the polarization vec- 
tor of the light field is oriented along the principal axis x 
or y, respectively. In the cases of parts (b), (c), and (d), 
where p — tt/6, 7t/4, and tt/3, respectively, the profile 
has double-peak structure. It is interesting to note that 
the two peaks have almost the same width. The reason 
is that, for the chosen parameters, the different princi- 
pal x and y polarizations of the light field correspond 
to almost the same reflectivity of the FBG mirrors and, 
consequently, to almost the same finesse of the cavity. 

We plot in Fig. 11 the tuning profile of the cavity 
transmission \T\ 2 for different values of the grating period 
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Fig. 11. Tuning profile of the cavity transmission |T| 2 = 
\T X \ 2 cos 2 ip+ \T y \ 2 sin 2 <^ for the grating period Ag = 363 
nm (a), 364.4 nm (b), and 366 nm (c). The orientation 
angle of the field polarization vector is (p = tt/4. The 
number of steps is N — 100. Other parameters are as for 
Fig. 8. 



Ag- It is clear from the figure that the tuning profile of 
\T\ 2 is very sensitive to the magnitude of the grating 
period Ag- In the case of part (a), where Ag = 363 nm, 
the mirror reflectivity |i?j,| 2 = 0.96 and cavity finesse 
F y = 70 for the y polarization are large as compared to 
the corresponding values l-R^ 2 = 0.16 and F x = 1.5 for 
the x polarization. This is the reason why a narrow peak 
appears on top of a broad peak. In the case of part (c) , 
where Ag = 366 nm, the mirror reflectivity and cavity 
finesse for the y polarization are negligible. Therefore, 
we observe in this case a significant background in the 
tuning profile of the cavity transmission \T\ 2 . 

5. Summary 

In conclusion, we have studied nanofibers with Bragg 
gratings from equidistant holes. We have calculated an- 
alytically and numerically the reflection and transmis- 
sion coefficients for a single grating and also for a cavity 
formed by two gratings. We have shown that the reflec- 
tion and transmission coefficients of the gratings substan- 
tially depend on the number of holes, the hole length, the 
hole depth, the grating period, and the light wavelength. 
We have found that the reflection and transmission coeffi- 
cients of the gratings depend on the orientation of the po- 
larization vector of light with respect to the holes. Such 
a dependence is a result of the fact that the cross section 
of the gratings is not cylindrically symmetric. Due to 
this property, the polarization vector of the field may be 
rotated after being reflected from or transmitted through 
the gratings. In addition, the tuning profile of the cav- 
ity transmission may have double-peak structure, whose 
components result from the resonances of the field with 
the different principal x and y polarizations. We empha- 
size that the only approximation used in our treatment 



is the omission of the radiation modes from the coupled- 
mode equations. This approximation is valid when the 
losses to the radiation modes are not serious, that is, 
when the characteristic size of the holes is small as com- 
pared to the fiber radius. It is interesting to note that, 
even when the length and depth of the holes are sub- 
stantial fractions of the fiber radius, our analytical and 
numerical results are in good qualitative agreement with 
the experimental results [29]. 
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Appendix A: Mode functions of the fundamental 
guided modes of a nanofiber 

For the fundamental guided modes, the propagation con- 
stant /3 is determined by the fiber eigenvalue equation 
[31] 



Jo(ha) 
haJiQia) 



n\ + n 2 K[ (qa) 
2n\ qaKi(qa) 



1 



2n 2 



n\k 2 



K'^qa) 
qaKi(qa) 

1 1 



h 2 a 2 

2 



q 2 a 2 



h 2 a 2 



1 V2 



(Al) 



Here the parameters h — (n 2 k 2 — /3 2 ) 1 / 2 and q — (/3 2 — 
ri^k 2 ) 1 / 2 characterize the fields inside and outside the 
fiber, respectively. The notations J n and K n stand for 
the Bessel functions of the first kind and the modified 
Bessel functions of the second kind, respectively. 

We present the mode functions of the electric and mag- 
netic parts of the fundamental guided mode fi = (u),f,l) 
in the forms £ M = e^e ilv and = h^e av , respec- 
tively [31]. The cylindrical components of the phase- 
shifted electric-part guided-mode function e^) are given, 
for r < a, by 



M 

"r 

M 



ic|^M[(i 



fC 



h J\{ha) 
2q Kijga) 
P Ji{ha) 



s)J (hr)-{l + s)J 2 (hr)}, 
s)J 2 (hr)}, 
Ji(hr), (A2) 



h Ji(ha) 
-K7i^H[(l--)Jb(ftr) + (l 



and, for r > a, by 

4*0 = iC[(l-s)K (qr) + (l + s)K 2 (qr)}, 



-IC[(1 - s)K Q (qr) - (1 + s)K 2 (qr)} 
fC^K^qr). 



(A3) 



Here the parameter s is defined as s — (l/q 2 a 2 + 
l/h 2 a 2 )/[Ji(ha)/haJi(ha) + K'^qa)/ qaKx{qa)]. Simi- 
larly, the cylindrical components of the phase-shifted 



10 



magnetic-part guided-mode function are given, for 
r < a, by 

+ (l + si)J 2 (/ir)]. 



/3ft Ji(fta) 
- (l + si)J 2 (ftr)], 

ujfio J\ (ha) 



(A4) 



and, for r > a, by 



= fic^i[(l - s 2 )K (qr) - (1 + a 2 )^2(gr)], 



^ = ifC ^i [{1 _ S2)Ko{qr) + (i + a2 )Jf 2 (gr)], 



/>) = l lC^sK 1 {qr). 



(A5) 



The power of the field in the fundamental guided mode 
fi is given by 

P^ = \J M^z ■ (e (M) x h<">*)] dr. (A6) 
It can be shown that 

x [J 2 (fta) + J^fta)] + (1 + s)(l + si) 
x [Jf(fta) - J 1 (ha)J 3 (ha)}} 

X [K 2 1 (qa)-K^qa)} + (l + s)(l + S2 ) 

x [K^K^-KHqa)]}. (A7) 

In the case where the coefficient C is real, we have the 
following symmetry relations: 

„(«,/,«) _ p {u-f,l) - „(«,/ -0 

» 



„(«,/,«) _ p (u),-/,o _ _ p («,/,-0 
e («>/.0 = _ e (« >-/.0 = e (^,/ -0 



and 
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